Probing Mott lobes via the AC Josephson effect 
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The alternating-current (AC) Josephson effect is studied in a system consisting of two weakly 
coupled Bose Hubbard models. In the framework of the mean field theory, Gross-Pitaevskii equations 
show that the amplitude of the Josephson current is proportional to the product of superfluid order 
parameters. In addition, the chemical potential-current relation for a small size system is obtained 
via the exact numerical computation. This allows us to propose a feasible experimental scheme to 
measure the Mott lobes of the quantum phase transition. 
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Introduction. Recent development of experiments al- 
lows the detection of the quantum phase transition in 
strongly correlated many-body systems [l[ . The relevant 
physics is captured by the Bose-Hubbard model, which 
describes the competition between the kinetic-energy and 
potential-energy effects. The Mott insulator (MI) to 
superfluid (SF) phase transition of the Bose-Hubbard 
model was described qualitatively using a mean field the- 
ory and realized in a gas of ultracold atoms [H . Gen- 
erally, for the MI to SF transition, the superfluid order 
parameter predicted by the mean field theory can not 
be measured directly. However, in this letter, we pro- 
pose a feasible realization to detect the phase transition 
and obtain the Mott lobes of the order parameter. The 
main idea relates to the well known alternating-current 
Josephson effect, which has been well investigated in the 
Bose-Einstein condensates (BEC) @, H, H, 0] ■ 

The Josephson junction is composed of two superfluid 
parts separated by a tunnelling barrier. When a constant 
chemical potential difference is applied, the AC Joseph- 
son effect occurs, correponding an oscillating particle flow 
through the barrier. If one part is fixed in surperfluid 
state, the phase of other part should determine the ampli- 
tude of the AC current. Inversely, it is possible to detect 
the state of the target part via measuring the Josephson 
current across the barrier. 

In this letter, through the well established mean field 
theory and spontaneous symmetry breaking mechanism 
for quantum phase transition, the analytical analysis re- 
veals that the amplitude of the Josephson current is pro- 
portional to the product of superfluid order parameters 
of the two parts. Through the measurement of the AC 
Josephson current, the Mott lobes can be obtained in- 
directly. Moreover, numerical simulations show that the 
profile of the Mott lobes for small size system is consis- 
tent with that in thermodynamic limit. This opens a 
possibility to detect the BEC in small quantum device. 
As an application we discuss an experimental realization 
in a cavity-atom hybrid system based on the recent work 
inRef. [|. 

Hamiltonian and the mean field theory. The setup is 
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FIG. 1: (Color online) Schematic plot of the lattice model 
used in this work. (A) is a Bose-Hubbard model with the on 
site repulsion U triggering the Mott transition; (B) is a free 
boson lattice model with a chemical potential fi which controls 
the density of bosons in lattice A; (C) is a thin contact surface 
between A and B which can be represented by a direct weak 
tunnelling across A and B or a free boson lattice model with 
an on site repulsion larger than u. 



depicted as the following Hamiltonian 
H = Ha + Hb + He , 
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Here cti and hi are the boson operators. In this pa- 
per, we consider lattices A and B having an identical 
structure with size N for simplicity. Each lattice site i 
(= 1, 2, N) corresponds to two positions in lattices A 
and B, respectively. Among them, i £ C lies on the con- 
tact surface. We consider the weakly coupled case g <C K. 
Fig. 1 is a schematic representation of the setup. 

We start our investigation from the mean field approx- 
imation. The order parameters can be defined as the 
expectation values of boson operators a\(t) and b\(t) in 
the Heisenberg picture by ignoring the fluctuations, i.e., 



X/J a>i (t) = (<Zi(t)) , ^ M (t) = (bi(t)) 



(3) 
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where the average is taken with respect to the ground 
state. In this paper, we use ai and bi to denote the boson 
operators in the Schrodinger picture. In order to study 
the time evolution of the order parameter, we make use of 
the Gross-Pitaevskii (GP) equations in the lattice [oL fioj] - 
The GP equations corresponding to the Hamiltonian (Q]) 
read 
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-k ^ i>aj + U \ip a .i\ 2 V'o.i - 9ii>b,u (4a) 

(iJ) 



^V'ba = -« ^ Ipb j + M^V ~ 9i^a,h (4b) 



where g\ = g for i G C and g\ — for i ^ C. 

When g is small enough, we can define solutions of Eqs. 
Ball and EH as 



(5) 



^a.i = (j>a (t) $ a ,i, A,i = ^b (t) <t>b,i, 



where $ a i and <E>bi are solutions for the ground states of 
Eqs. (f4"a[) and (|4b|) when g = 0. $ i and are se t to 
be real and normalized as <1> 2 ; = i ~ 1- Taking 

the periodic boundary condition, $ a i = <&£,i = l/\/A. 
After replacing ^ 0j j and V^i, GP equations Ba)l and (|4b|) 
become the two-state model 111! 
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for two d-dimensional lattices with the periodic boundary 
condition and contact surface with size N~a~ . (j) ab (t) = 
\jN a ^e Ma - b where N a .b and 9 a ,b are the particle number 
and phase of the superfluid component in lattices A and 
B. 

Substituting tp a ^ (t) into the above two-state model, 
we get 



^ = -2A^1-^ 2 sine, 
at 



80 
dt 

where 9 = 9 a 



AE + Az - 
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%, z = (N b - N a )/(N a + N b ), AE 



E b -E a -\ (N a + N b ), and A = U a (N a + N b ) /2. Eqs. 



(16a)) and (|6b[) show that N a + N b is conserved. Then the 
Josephson current (if the boson is neutral, the current 
corresponds to the flow of bosons) is 



J(t) 



dt 



-2Ky/N a N b smQ. 



When fj, 3> U, k, it becomes 
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sin fit. 



(9) 



(10) 



Here N aib is time-dependent. On the other hand, the up- 
per bound of the particle immigration across the junc- 
tion during a half period of oscillation, which corre- 
sponds to the case N a , b ~ A, is of the order AJVj ~ 
g^flVJV b ~/(pN : 3) ~ gN^r/fi. For A^ < A, A a , fc can 
be regarded as a constant, which corresponds to the or- 
der parameter defined in the framework of the mean field 

theory | as </> 7 = ( 7i ) = (7/) = y/EJN, (7 = a,b). 
Then the Josephson current is obtained as 



J(t) = — 2gN d ipai'b sin (j,t. 



(11) 



Measurement of the order parameter. In this section, 
we focus on the experimental scheme to obtain the Mott 
lobes of the Bose-Hubbard model. It is based on mea- 
suring the magnitude and frequency component of the 
Josephson current. Our proposal for the measurement 
of the order parameter experimentally is as follows: (a) 
Parameters U, k and \x are set to be values correspond- 
ing to a point (n/U , k/U) on the phase diagram, then 
cool the system to the ground state, (b) Take the ground 
state as an initial state, and shift the chemical potential 
to fj, + A/z (A/i ^> U,k). (c) Obtain the current J(t) be- 
tween A and B. (d) Numerically analyze the curve J(t) 
to obtain 

J m (n/U, k/U) = max {J (ji/U, k/U, t)} (12) 

and the frequency component of the current from the 
formula 



J((jl/U,k/U,(jj) = \ - dtJ(n/U,K/U,t)smuit (13) 



as t — > 00. (e) In the weak coupling limit g -C k, all the 
bosons in lattice B are in the superfluid state, i.e., ip b ~ 
(blh) 1 / 2 . Then can be measured via the measurement 
of the average particle density in lattice B. (f) Construct 
the phase diagram ip a (fi/U, k/U) according to Eq. (TTTj) . 

Theoretically, the Josephson current arises from the 
time evolution of the superfluid state of the Hamilto- 
nian (fTJ) . In this Hamiltonian, we add the chemical 
potential /i only in system B, but this is equivalent 
to adding a chemical potential —fx in system A, since 

adding a term -fJ,J2i ( a i°i + ^t^i) to the whole system 
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FIG. 2: (Color online) 3D plots of the Josephson current (a) 
J (U,t) and its Fourier component (b) J (U,ui) for a small 
size system obtained by the exact diagonalization with k = 1, 
Afi = 100, g — 0.1 and n/U — 0.5. In (b), the upper limit 
of the time integration is taken as t(1/A/i) = 20. As r in- 
creases, the peak with frequency 100 will become higher and 
sharper. The current drops rapidly as U increases. These 
results show that the current is a Josephson alternating cur- 
rent with frequency Afi and can be employed to witness the 
quantum phase transition. 

does not bring any physical change. Then the order 
parameter of lattice A corresponds to the point (fJ./U, 
k/U) of the obtained phase diagram. At t = 0, the 
ground state of the Hamiltonian {1} is \<p g (n/U, k/U)) = 
\tp(0)), which is the initial state for the time evolution 
driven by the Hamiltonian H + A/i £\ b\bi, i.e., \(p(t)) = 
g-ifH+ApE,^,) \ip(0)). Then the current across A and 
B is 

J(li/U,K/U,t) (14) 
= -zgJ2(v(0)\(4(t)bi(t)-h.c)\<p(0)), 

where a\ (t) and bi (t) are boson operators in the Heisen- 
berg picture. To demonstrate this scheme, the numerical 
simulation is performed for N = 2 system with k = 1, 
A/i = 100, g = 0.1, and fi/U — 0.5, which corresponds to 
a parallel line in the phase diagram of lattice A. Fig. 2(a) 
is the 3D plot of the Josephson current J (U, t), which is a 
good sinusoidal curve vs time t when U is not very large, 
in which region systems A and B are both in the super- 
fluid phase. When U is large enough, the amplitude of 
the current drops, which indicates that system A enters 
into the Mott insulating phase. Fig. 2(b) is the 3D plot 
of the Fourier component J (U, to) obtained according to 
Eq. (Tl3"|) . The upper limit of the time integration (fl"3"|) 
is taken as t(1/A/x) = 20. The larger values of r will 
make the peak with frequency 100 higher and sharper. 
The current drops rapidly as U increases, which indicates 
the quantum phase transition of system A. These results 
show that the current is a Josephson alternating current 
with frequency A/i and can be employed to witness the 
quantum phase transition. 

According to the mean field results, when g is small 
enough, the Josephson current is proportional to the or- 
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FIG. 3: (Color online) Mott lobes calculated by the exact di- 
agonalization with an auxiliary field for a small size system 
(color contour map), and by the mean field method via the 
Gutzwiller trial wave function. The contour line denotes the 
boundary of two phases which corresponds to a vanishing or- 
der parameter. It is shown that although the result of the small 
system can not give the boundary of two phases, its contour 
lines are consistent with that from the mean field method in 
the thermodynamic limit very well. 

der parameter of an isolated system with the Hamilto- 
nian Ha — > Ha — MX)i a i a i- However, the Mott lobes 
obtained in a finite size system from the above proce- 
dures (a-f) cannot be compared with the phase diagram 
obtained from the mean field method in the thermody- 
namic limit, which is usually obtained via the Gutzwiller 
trial wave function 0] • To compare the Josephson current 
with the mean field phase diagram, we need to develop 
another way to obtain the order parameter for a finite 
system in the framework of the spontaneous symmetry 
breaking. It is well known that the quantum phase tran- 
sition is a consequence of U(l) symmetry breaking for 
an interacting boson system [12j. Then introducing an 
auxiliary field as 

HA~li^2 a l ai ~ * Ha—(*'5~* j a^ai+A^ (a\ + aij (15) 

i i i 

breaks the U(l) symmetry breaking and induces the or- 
der parameter. Accordingly, the order parameter is de- 
termined by ip a = lim,\^o limjv-«x> ( a i) , where the av- 
erage is taken for the ground state of the Hamiltonian 
(|15p . In the thermodynamic limit N — > oo, the ground 
state \<Pg) is a coherent state and the vanishing ip a dis- 
criminates two phases. In this case, the order parameter 
is independent of A. On the other hand, when lattice A 
is coupled to B, the order parameter obtained from Eqs. 
(TrTj) and (TT4")) should also be independent of g in the limit 
N — > oo and g — > 0. 

In a finite system, although ijj a obtained from the above 
two ways are not independent of A and g, it is believed 
that their consistency should be revealed from the con- 
tour maps. To demonstrate this, the numerical simula- 
tion for a small size system is performed and compared 
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FIG. 4: ( Color online) Mott lobes of the order parameter for a 
small size system calculated by the exact diagonalization with 
K = 1, A/i = 100, and g = 0.1. The color contour map is 
obtained from the Josephson current, while the contour lines 
are obtained via an auxiliary field with A = 0.1. It is shown 
that two results are consistent very well. 



with the mean field method. Even for a small size sys- 
tem, the dimension of the Hilbert space of the Hamil- 
tonian (|15[) is infinite. So the truncation approximation 
is taken for the exact diagonalization of the matrix for 
N = 2 and the average density (a\aij £ [0,2]. In Fig. 

3, the Mott lobes calculated by the exact diagonalization 
with an auxiliary field for N — 2 system and by the mean 
field method via the Gutzwiller trial wave function are 
plotted. The contour line denotes the boundary of two 
phases, which corresponds to a vanishing order param- 
eter. It is shown that although the result of the small 
size system can not give the boundary of two phases, its 
contour lines are consistent with that of the mean field 
method in the thermodynamic limit very well. 

In Fig. 4, Mott lobes of the order parameter for a 
small size system arc plotted through two different mech- 
anisms with k = 1 and g = 0.1. The color contour map is 
obtained from the Josephson current via the above pro- 
cedures (a-f), while the contour lines are obtained from 
the Hamiltonian (| 1 5|) by the exact diagonalization. It is 
shown that two results are consistent very well, which 
indicates that the property of the small size system can 
shed light on the profile of Mott lobes in the thermody- 
namics limit. 

Discussion. In order to detect the Mott lobes of the 
quantum phase transition, we investigate the AC Joseph- 
son effect in a system consisting of two weakly coupled 
Bose-Hubbard models. The mean field theory in the 
thermodynamic limit and the numerical simulation for 
a small size system show that, through measuring the 
magnitude and frequency component of the Josephson 
current, the Mott lobes can be measured. 

To realize this scheme experimentally, as mentioned 



before, a good candidate is the coupled cavity system 
with each cavity interacting with 4-level atoms driven by 
an external laser Q. In this system, the repulsion U 
can indeed be strong enough to observe the Mott insu- 
lator state for photons. Moreover, according to Refs. [§] 
and [l3| , the controllable range of the chemical potential 
required by our scheme is in experimentally accessible 
parameter regimes. In fact, in such an effective Bosc- 
Hubbard system, the on-site interaction U and chemical 
potential /i for photons are determined by 



u = s 



(913\ 2 924 

V n J a : 



(16) 



under the conditions g, A, 524, e <C ^; <?24ff |Af2|. Here 
f2 is the Rabi frequency of the driving laser; S is the atom 
number in each cavity; 913 , 924 are the couplings between 
cavity mode to the atomic levels; S, A and e are detunings 
of atomic transitions with respect to the cavity and laser 
fields. All the notations are identical with those used in 
Fig. 1 of Ref. 0. 

By fixing S (gi 3 /Q) 2 and adjusting g| 4 /A and e, U and 
\x are tunable independently. This allows us to simulate 
the Hamiltonian ([2]) by setting e = in lattice A to get 
li = 0, and (724 = in lattice B to get the vanishing U . 
Subsequently, to drive the quantum phase transition and 
probe the Mott lobes of lattice A, the parameters in two 
systems A and B are tuned according to the procedures 
(a-f). In this scheme, the Josephson current J(t) is the 
flow of photons in photonic crystal waveguides, which 
can be imaged via a high-resolution imaging technique, 
the collection scanning near-field optical microscope [l4( . 
This predicts that the Mott lobes for a small quantum 
device can be detected experimentally. 
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